Generating non-Gaussian states using collisions between Rydberg polaritons 



(N 

o 

(N 



O 
m 



Oh. 



> 

rn 

O 



X 



Jovica Stanojcvic"'^, Valcntina Parigi^. Erwan Bimbard^. Alcxci Ourjoumtsev^ , Pierre Fillet^ and Philippe Grangier^^ 

^ ^ Laboratoire Charles Fabry, Institut d'Optique, CNRS, 

Univ Paris-Sud, 2 Avenue Fresnel, 91127 Palaiseau, France. 

^ Laboratoire Aime Cotton, Batiment 505, Univ Paris-Sud, 91405 Orsay cedex, France 

(Dated: March 1, 2013) 

We investigate theoretically the deterministic generation of quantum states with negative Wigner 
functions, by using giant non-linearities due to coUisional interactions between Rydberg polaritons. 
The state resulting from the polariton interactions may be transferred with high fidelity into a 
photonic state, which can be analyzed using homodyne detection followed by quantum tomography. 
Besides generating highly non-classical states of the light, this method can also provide a very 
sensitive probe for the physics of the collisions involving Rydberg states. 



The generation and characterization of highly non- 
classical states of the light have accomplished consider- 
able progress during recent years. This includes for in- 
stance the production and analysis of one- [l[ and two- 
0, 9 photon Fock states, and of superpositions of co- 
herent states, often called "Schrodinger's cat" states [J]. 
In these experiments, the desired states are obtained by 
using so-called "measurement-induced non-linearities" , 
where a measurement is performed on one part of an 
entangled state. Then the other part is projected onto 
the desired state, conditional to obtaining the good re- 
sult in the projecting measurement. This method works 
quite well, but it is intrinsically non-deterministic : the 
probability of success is usually low, and the desired state 
cannot be created "on demand" , when needed for appli- 
cations e.g. for quantum communications. 

Here we would like to investigate another scheme, 
which can, at least in principle, be made more deter- 
ministic, by using Rydberg interactions in a cold atomic 
gas. Rydberg states are highly excited atomic states, 
that interact very strongly at distances R of order a few 
/xm, through either dipole-dipole {l/R^) or van der Waals 



{1/R ) interactions 



ipoie- 



The idea is first to change 



a generic photonic state, for instance a weak coherent 
state, into a so-called "polariton" state, where a long- 
lived excitation is distributed over many atoms [7[. Here 
we will consider Rydberg polariton states, where many 
atoms share a few delocalized Rydberg excitations Q- 
The prepared polariton state will then evolve under the 
action of these Rydbcrg-Rydberg interactions, generat- 
ing highly non-classical (typically non-Gaussian) states 
Q. When the desired state is obtained, the polariton 
state can be converted back into a free-propagating pho- 
tonic state, using a control laser beam (see Fig. [T|). The 
phase-matching condition between the input, write, read, 
and signal beams leads to a collective enhancement effect 
ensuring that the state of the light is emitted in a well- 
defined spatial and temporal mode [8| . The conversion of 
the polariton state back to photons is then expected to 
deterministically produce highly non classical states |10l - 
Il3j . In addition, the analysis of the generated state will 
give information about the coUisional mechanisms which 
take place between the Rydberg atoms. 
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FIG. 1: An ensemble of A'^ three-level atoms is initially pre- 
pared in a coherent Rydberg polariton state, with only a few 
atoms excited. After a time shorter than the decoherence 
time, the collisions have modified in a non-linear way the co- 
efficients of the polariton state, which is "remapped" onto a 
photonic state, and analyzed using an homodyne detection. 



Let us emphasize that a particularly useful method to 
fully characterize highly non-classical states is quantum 
tomography, which allows one to reconstruct the Wigner 
function W{q,p) of the state in phase space, where q and 
p are the quadrature operators of the quantized electric 
field, measured using homodyne detection [ll-Q]. Such a 
method provides a full characterization of the measured 
quantum state, which is very intuitive because the "non- 
classicality" of the measured state, and especially its pu- 
rity, directly translate into the negativity of the Wigner 
function [ij] . We will thus use such Wigner functions to 
characterize both the photon and the polariton states. 

In this letter we will study the first two steps of such 
a deterministic preparation of non-Gaussian states : 

(i) preparing the coherent polariton state using a weak 
laser pulse; we will introduce the phase-matched symmet- 
ric Dicke states [lH as a convenient way to describe the 
states of the ground and excited atoms (this is essentially 
equivalent to the polariton picture [7|). 

(ii) leaving the polariton state evolve for some time 
(shorter than its decoherence time). The state will then 
be modified in a non-linear way, depending on the na- 
ture of the Rydberg-Rydberg interactions. We will ob- 
tain simple analytic expressions for this evolution, that 
are the main results of the present paper. 

Finally, we will characterize the generated polariton 
states by computing their Wigner functions, and we will 



discuss various experimental considerations, including 
the remapping of the polariton into a photonic state. 

The calculation is performed by splitting the evolution 
of the system in two steps : first, an excitation step using 
a short (typically ~ 1 ns) and weak laser excitation pulse, 
creating a few Rydberg states; second, a free evolution of 
the generated state under the effect of Rydberg- Rydberg 
interactions. Wc will show that we can consistently ig- 
nore the interactions during the first step, and then ig- 
nore the laser in the second step, since it is turned off. 

In order to describe the excitation step for an ensem- 
ble of N atoms from a ground state g, we consider a 
two-photon excitation, off-resonant from the intermedi- 
ate level r, and resonant with the Rydberg state e (see 
Fig. 1). It can thus be described using an effective two- 
level model m with the Hamiltonian H = Hi+Hc, where 
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Here ctJ,^ = \a) (^| e±(i-'5°'3)*'^-R', R, is the position of 

atom i, and k the total wave- vector of the exciting light, 
with an effective (pulsed) Rabi frequency ri(i); fiKij is the 
pair- wise interaction energy of the excited atoms, and 5a/5 
the Kroneckcr symbol. 

For short times and very low excitation fractions, let 
us first neglect the interaction term He and keep only 
the laser excitation Hg. The result of laser excitation is 
then a coherent polariton state, that is a superposition 
of symmetric "phase-matched" Dicke states \n), where n 
is the number of excited atoms (see Appendix 1). The 
amplitude C„ corresponding to the collective states \n) 
directly follows from the single-atom amplitudes 



Cn = (-i)"v/5(A^,n)sin" - cos 



N-n 



U! 



(1) 



where B{N,n) is the binomial coefficient, A^ is the total 
number of atoms and w = J n{t)dt is the pulse area. In 
the limit of very small n/N, one has 



lim Cn/C„ 



(2) 



where a = — iviVtan(aj/2) w ~iy'Nuj/2. The amplitude 
a is related to the averaged number of excited atoms 
(iVoxc) = Nuj^/4: = |a|2. Relations © yields the well 
known expression Cu/Cq = a"/vn! for the amplitudes of 
a coherent state, with Cq determined by normalization. 
For this calculation to be consistent, we need to fulfill 
the condition {n\ exp {—iHcT/h)\n) ^ 1 at the end of the 
laser pulse of duration r, where He is the collisional part 
of the Hamiltonian; this will be checked in Appendix 2. 
After the laser pulse is off, we consider the evolution 
of the previous coherent Rydberg polariton state under 
the only action of the Rydberg interaction Hamiltonian 
He- A crucial remark is that during the excitation and 
de-excitation phases (see Fig. 1), only the the symmet- 
ric phase- matched Dicke states (see Appendix 1) will be 



mapped coherently between the polariton and photonic 
states [7| . In addition, the Hamiltonian H^ preserves the 
number of excitations n, and therefore it only mixes sym- 
metric Dicke states with non-symmetric ones, which are 
uncoupled from the laser readout process. 

In order to characterize the phase-matched part of the 
state after an evolution time T, wc need therefore to eval- 
uate the matrix elements {n\U\n), where U = e~^^''^/'^. 
For this purpose wc use the following transformation 



= 1 + a^'' ct'' 



-1 



(3) 



For a low number of excited atoms, the probability to 
have p interacting atoms close to each other vanishes very 
quickly as p increases. Therefore, the leading interaction 
order originates from two-body interactions, the next- 
to-leading order originates from three-body interactions 
and so forth. The transformation ([3]) can facilitate this 
expansion because the term (e~*''p9^ — 1) is zero if the 
two atoms do not interact, i.e. if they are far from each 
other. One can then use these types of terms to select 
pairs of interacting atoms in various interaction orders, 
and the expectation values {n\U\n) can be evaluated by 
bookkeeping various combinations and contributions of p 
excited atoms, with p < n. The expansion is finite (since 
p<n), however the number of terms and their complexity 
rapidly increase for p > 4; we will therefore look first at 
low n, and then find an excellent ansatz for higher n. 

Denoting as ri[r) the atom number density at point r, 
and dNi ~ d^ri i]{vi), we define : 



J3 == ^ / / / dN, dNi dNs X 
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(1|C/|1) = 1, 
I + /2, 

I + 3/2 + 3J3+ /3, 

l + B(n,2)J2 + 3S(n,3)J3 



(e-'«..'r_l)(e-"^"^-l)(e- 

and one gets the successive terms 

(0|C/|0) 
(2|C/|2) 
(3|C/|3) 
{n\U\n) = l + B(n,2)J2 + 3S(n,3)J3+ B{n,'S)l3 + ■ ■ ■ 

For each {n\U\n), the coefficients of the quantities /2,3 
and J3 correspond to the numbers of different choices of 
excited atoms {i,j) and (i,j, s) that appear in the ex- 
pressions (H])-®. Though this approach is rigorous and 
can work in principle at any order, it has the disadvan- 
tage that the integrals /„ are more and more complicated 
to evaluate when n increases. Whereas I2 can easily be 
calculated analytically (see Appendix 3), this is more te- 
dious for J3, and /„>3 are only numerical. We therefore 
introduce now a much simpler approach, giving analyt- 
ical results at any order, which works surprisingly well 
when compared with numerical calculations. 



The idea is to evaluate Un>m , assuming that Um ~ 
{m\U\m) is known. For this purpose, we first note that 
the value of [/„ for a set {n\ = {ii, ^2, ... i„} of n 
excited atoms can be formally divided into a product of 
terms involving set {m\ oim < n atoms, in the following 
way (we note that the total number of excited atoms is 
always conserved) : 



{cxp{-iKc))[,-,y= Y]_ (<2xp(-JA„mA', 



(7) 



{'"} 



where K^ = HcT/h, and X„rH^B{n,2)/[B{n,m)B{m,2)] 
is due to the fact that each k^ appears in several subsets. 
We also note that the calculation of Um for a given set 
{n} of atoms will involve an average over all random po- 
sitions of atoms, that is essentially the continuum limit 
of the above expressions. Taking into account this aver- 
aging, the quantities (exp(— iA„m-f'^c)){m} are the same 
for all r7i-subscts. As a consequence, we use as a crucial 
ansatz that the rhs of Eq. ([7|) is just a product oi B{n, m) 
identical factors (cxp(— iA„m-?''c)){m} • 



{n\U\n) = {exp{-iKc))^„y 

« (exp(-zi^,))f^l"j-")'"'" = (m|C/|m)^^^fe^ 

For 771 = 2, one has then simply: 

(n|[/|n) = (l + /2(r))"("-i)/' 



(8) 



(9) 



where I2 (T) can be calculated analytically by integrating 
on the positions of two atoms within a sphere of constant 
atomic density, the result is given in Appendix 3. 

Figure [2] shows the plots of (n|C/|n), n = 2, 3 and 4 for 
van der Waals interactions, as a function of the scaled 
interaction time t. For each n three different curves are 
almost perfectly overlapping : the green dotted curves 
are the numerical results, the red full ones show the scal- 
ing from n = 2 using the expression of l2{t) given in 
Appendix 3, and the blue dashed ones show the scaling 
from n = 3. In the numerical calculation, groups of four 
atoms are generated and {n\U\n) is essentially the aver- 
aged U over random n-groups of atoms. The calculations 
are done for a sphere with a uniform density, but can be 
generalized for arbitrary density profiles. 

Summarizing, we have thus obtained a series of sim- 
ple approximate expressions of {n\U\n), valid for any n. 
This surprisingly simple derivation can be understood as 
an approximate but efficient way to resuni the terms ap- 
pearing in the more rigorous expansion quoted before. 

In order to analyze the generated polariton states, it is 
convenient to use Wigner functions, that show the evo- 
lution of the initial Gaussian into non-Gaussian states. 
Simple analytical calculations give easily W{q,p) in a 
suitably truncated Fock state basis, as a function of the 
coefBcients {n\p\m) of the density matrix, obtained from 
the evolution of the initial Dicke state using previous for- 
mulas. For long interaction times, the result is rather 
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FIG. 2: Real and imaginary parts, and modulus of Um ~ 
{m\U\m) as a function of time in units of Tr — _R^/|C6|, for 
attractive van der Waals interactions. For repulsive interac- 
tions one has only to change the sign of Im(C/„). All red curves 
are obtained by scaling from the corresponding U2 using Eq. 
®. All blue dashed lines are obtained using the scaling for- 
mula ((8| with the exact numerical result for C/3 , and all green 
dotted lines are the results of the full numerical calculation. 
Similar results (not shown here) can be obtained for dipole- 
dipole interactions rather than Van der Waals interactions. 



simple : since the coefficients for n = and 1 are unaf- 
fected, and all other ones go to zero, the Rydberg medium 
behaves as deterministic near-perfect "quantum scissors" 
[l8|, cutting the initial coherent state into the subspace 
of Fock states with zero and one photon : 



a 



\0)+a\l) 



~^ For the (repulsive) 70s state of Rb*^, 
= 880 GHz fimP [l3| , so the time scal- 



This is in agreement with the recent result of Kuzmich 
et al. [9| . In addition, we have the explicit expressions to 
calculate the evolution between the initial Gaussian co- 
herent state, and the final truncated non-Gaussian state, 
as illustrated on Fig. [3] 

As an cxemplc of realistic experimental parameters, let 
us consider 2500 atoms in a sphere of radius i? = 10 /im so 
Uat = 6 10^^ cm" 
one has \Ce\/{2TT) 

ing is Tfj = R^/\Ce\ = 180 ns. An initial coherent state 
with amplitude |a| ^ 1, excited by a pulse of duration 
0.72 ns (0.004 Tr) will evolve into a non-Gaussian state 
within a time of 720 ns (4Tr), as shown on Fig. [3l If 
the atoms are inside a low-finesse optical cavity, these 
numbers also warrant a large cooperativity parameter C 
with high outcoupling mirror transmission Tc (C ~ 180 
for Tc ~ 0.01), which should in turn warrant a good re- 
covery of the photonic state [10, [ill, [i3| ■ 

As a conclusion, we have studied the evolution of a co- 
herent (Dicke) Rydberg polariton state, under the effect 
of Rydberg- Rydberg collisions. The non-linearities are 
clearly large enough to have an effect at the few-photon 
level, even outside the dipole blockade range. They are 
able to turn an input coherent state into a non-Gaussian 




FIG. 3: Evolution of the Wigner function of an initial coher- 
ent state \a) with \a\ = 1, between the scaled times 0.004 (left 
side) and 4 (right side), that is close to the truncated state 
|0) +a|l). The corresponding physical times are 0.72 ns (left) 
and 720 ns (right), with Van der Waals interactions only. 



state, at the expense of significant losses due to the de- 
coherence of states containing more that one polariton. 
Whether or not such decoherence effects can be avoided 
in order to reach a high input-output recovery of photonic 
states is still an open question. Let us emphasize how- 
ever that this decoherence docs not prevent the polariton 
state, once created, to be remapped with high efficiency 
on a photonic state, and then analyzed using an homo- 
dyne detection [13| . This may also provide an interesting 
way to investigate Rydberg-Rydberg collisions. 



Appendix 

(Al) The Dickc state \n) corresponding to n ex- 
cited atoms among TV atoms is defined as the eigenstate 
with eigenvalues (n — N/2, r{r + 1)) of the operators 



?(^) 



(J^^^ j|), with Ji 

?(+) _ v-W 
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^N 



= E»=i(|e,*>(e,i| - |.g,i)(5,i|)/2, 



= e; 



e,i)(5,i|exp(ifc.i?i), J\ 



(-) _ 



?(.+)^t 



.j(-)j(+) 



{JY>y and 
)/2 + (jl^^)2, where R, is the 



k ^ ^^i=l 

position of atom i, and k is the total wave- vector of the 
exciting light. Symmetric Dicke states are obtained for 
r = N/2, and non-symmetric ones for r < N/2 — 1. 

(A2) To estimate the action of the collisions during the 
short excitation phase of duration t, two points must be 
considered. The first one is that the excitation of a Ryd- 
berg atom creates a "blockade sphere" around it, where 



a second atom cannot be excited. The second one is that 
the quantities U„ calculated above should remain close to 
one during r, so that the excitation and the collisions act 
on different time scales. To evaluate the first correction 
we can use the following correlation function Gij derived 
in [131 for two level atoms and small a; 



Gij = 



J[ dhe'^^^^^n{ti)u}{ti) 



.){^i 



Hrr 



(10) 



This function prevents two excited atoms to be closer to 
each other than a distance r^ with magnitude given by 
Cgr/r^ ^ 1, where r is the duration of the excitation 
pulse. It is thus clear that r^ will be small if r is small. 
More quantitatively, the modification of Eq. (jH]) would 
be to substitute /2(r) by c{t,T)I2{T), where c{t,T) is a 
correction factor evaluated numerically. For the parame- 
ters quoted in the text, we numerically get c(t,t)~ 0.984 
at the end of the excitation pulse, and this effect is thus 
negligible. Similarly, using the results in Appendix 3 one 
gets for t = t/Tr = 0.004 the values U2 w 0.94 + 0.05 i, 
J/a w 0.82 -H 0.13 i, Ui « 0.66 + 0.21 i, that is acceptably 
close to one if \a\ is small enough (see also Fig. [5|). 

(A3) The analytical expression of I-iit) used in the 
curves on the previous page is : 



hit) 




it 
3'~64 



where the second expression is a short-time approxima- 
tion valid up to i ~ 0.1. Here f is a scaled time related to 
the physical time T by i = T/Tr, where Tr^ R^/\C(i\, 
and R is the radius of the spherical volume of the sample. 
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